Distortion of the Stoner-Wohlfarth astroid by a spin-polarized current 
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The Stoner-Wohlfarth astroid is a fundamental object in magnetism. It separates regions of the 
magnetic field space with two stable magnetization equilibria from those with only one stable equilib- 
rium and it characterizes the magnetization reversal of nano-magnets induced by applied magnetic 
fields. On the other hand, it was recently demonstrated that transfer of spin angular momentum 
from a spin-polarized current provides an alternative way of switching the magnetization. Here, 
we examine the astroid of a nano-magnet with uniaxial magnetic anisotropy under the combined 
infiuence of applied fields and spin-transfer torques. We find that spin-transfer is most efficient 
at modifying the astroid when the external field is applied along the easy-axis of magnetization. 
On departing from this situation, a threshold current appears below which spin-transfer becomes 
ineffective yielding a current-induced dip in the astroid along the easy-axis direction. An extension 
of the Stoner-Wohlfarth model is outlined which accounts for this phenomenon. 
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I. INTRODUCTION 



The study of magnetization reversal has provoked con- 
tinuous interest during the last sixty years. Indeed, this 
topic has proven to be scientifically very challenging, and 
one of the fundamental issues in magnetic data stor- 
age and memory technologies. The difficulty to quan- 
titatively understand magnetization reversal, known as 
Brown's paradox,^ results from complex domain patterns 
that can form in magnetic materials. However, for nano- 
scale systems the picture is greatly simplified, as ini- 
tially described by Stoner and Wohlfarth^ and Neel?^ 
where the magnetic order parameter is assumed to be 
uniform across the dimensions of the sample and can be 
described as a single macro-spin. Fifty years were nec- 
essary to have the technological environment needed to 
test the macro-spin model experimentally During that 
time, new techniques essential to fabricate and character- 
ized nanometer-size objects such as scanning probe mi- 
croscopies and nano-lithography were developed. It was 
then demonstrated that, for some magnetic nanopartic- 
ules, the angular dependence of the switching field can 
indeed follow the famous Stoner-Wohlfarth astroid curve 
which separates regions of the cartesian magnetic field 
space with two stable magnetic states from those with 
only one stable state. 

Very few magnetic systems, however, possess the pure 
uniaxial magnetic anisotropy required to closely obey 
the Stoner-Wohlfarth model and the locus in the field 
space of all the magnetization switching fields of small 
magnets seldom take the shape of an astroid, as defined 



mathematically. However, over the years, the concept of 
Stoner-Wohlfarth astroid has been generalized to arbi- 
trary effective anisotropyi^ The term is now commonly 
used to refer to any critical surface that delimits the 
region of bistability of the magnetization in the field 
space and is a fundamental property of magnetic nano- 
materials. 

More recently, it has been theoretically predicted— 
and experimentally evidenced— that injection of spin- 
polarized electrical currents can induce magnetization 
reversal. This new approach for magnetization rever- 
sal has generated considerable scientific interest and is 
expected to play a major role in many emerging spin- 
tronic technologiesi^ Usually, transfer of angular momen- 
tum from a spin-polarized current to the magnetization 
does not strongly modify the positions of the magneti- 
zation equilibriaF^^ Its primary effect is to produce an 
extrinsic damping which either reinforces or opposes the 
intrinsic damping of the magnetization and modifies the 
stability of the equilibria rendering, for example, an un- 
stable equilibrium stable. Spin transfer is expected to 
produce significant distortions of the astroid, making the 
region of bistability expand in some parts of the magnetic 
field space and retract in others. This distortion was con- 
sidered theoretically by Sun— in the case of an in-plane 
polarized current acting on a magnet with a combination 
of uniaxial and easy-plane magnetic anisotropics. 

In this paper, we examine how the Stoner-Wohlfarth 
astroid is indeed distorted by spin-transfer in the simplest 
and most fundamental case of a nano-magnet with uniax- 
ial magnetic anisotropy submitted to a current polarized 
along the easy-axis. As we do this, we reveal how the ef- 
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ficiency of the spin-transfer torque with respect to mag- 
netization switching varies in a rather counter-intuitive 
manner with the orientation of the external magnetic 
field. 



II. EXPERIMENTS 

The small magnet whose astroid is studied experimen- 
tally is the 3-nm thick 'free' element of a giant magne- 
toresistance (GMR) spin-valve patterned in the shape of 
a vertical pillar with 100 nm x 200 nm hexagonal cross- 
sectional shape. Both the free-element and the reference- 
element of the device have perpendicular-to-plane mag- 
netic anisotropy. This is achieved by using [Co/Ni] and 
[Co/Pt]/[Co/Ni] multilayer stacks, respectively, as de- 
scribed in details elsewhere.— In comparison with the in- 
trinsic perpendicular-to-plane anisotropy of the extended 
Co/Ni multilayer, the shape anisotropy introduced in the 
plane by patterning the film in a non-circular pillar is 
negligible. Thus, the overall magnetic anisotropy of the 
free-element is essentially uniaxial and its astroid should 
be in first approximation invariant upon rotation around 
the normal to the film plane z. 

In practice, the astroid was built by determining the 
switching fields of the free-element from differential re- 
sistance dV/dl versus magnetic field minor hysteresis 
loops [Fig. [Ha)] recorded at varying field angle Oh from 
z (-90° <eH < +90°), in the plane defined by z and the 
small axis of the hexagon y (see inset in Fig. [Ijb)). The 
differential resistance was measured using a Lakeshore 
Model 370 AC resistance bridge with an excitation cur- 
rent of 10 AiA RMS at 13.7 Hz. The DC current provided 
by a Keithley Model 2400 sourcemeter was injected in 
the sample using a home-made bias-T interface. Before 
measuring each magnetoresistance loop, a positive field 
in excess of 1 T was applied to ensure that both the 
free-element and the reference- element would be initially 
magnetized positively, i.e. in the z > semi-space. The 
experiments have been performed several times and no 
strong fluctuation of the switching fields, i.e. a stochastic 
behavior related to thermal fluctuations, was observed. 
Moreover, no signature of spin-transfer-induced steady 
precession states, i.e. non-hysteretic peaks or dips in the 
differential resistance versus field loops^^ was found with 
the DC current values used. 

In the off-axis geometries {Oh 0), the magnetiza- 
tion vectors of the two elements do not always remain 
strictly colinear to each other during field cycling. This 
manifests in the curvature of both the decreasing-field 
branch and, more obviously, the increasing- field branch 
of the measured GMR loops [Fig.[IJa)]. However, for the 
sake of simplicity, we will still name the low-resistance 
state and the high-resistance state of the nanopillar the 
'parallel' (P) state and the 'antiparallel' (AP) state, re- 
spectively. Three dimensional micromagnetic simula- 
tions performed with the OOMMF software packageJ^ 
indicate that the continuous change in the relative ori- 
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FIG. 1: (a) Examples of differential resistance versus mag- 
netic field minor loops of the free element in a 100 nm x 
200 nm [Co/Ni] /Cu/ [Co/Ni] /[Co/Pt] nanopillar recorded at 
diff'erent angles of the applied field with respect to the film 
normal. For clarity, the curves corresponding to field an- 
gles other than 0° are offset vertically, (b) Experimentally 
determined two-dimensional cross-section of the zero-current 
astroid of the free-element, in the {Hy,H^) plane. Circles 
indicate a transition from the parallel (low-resistance) state 
to the antiparallel (high-resistance) state of the GMR device. 
Squares indicate the inverse transition. 



entation of the magnetization vectors that the curvature 
reveals is primarily ascribable to the rotation of the mag- 
netization of the free-layer. To a good approximation, 
the magnetization of the reference-layer, harder magnet- 
ically, remains fixed. An unambiguous determination of 
the switching fields Hap^p and Hp^ap was possible 
only for field angles such that \0h\ < +85°, where the re- 
versals of the free-element magnetization occur abruptly. 
For 85° <\0h\< 88°, the GMR loops may contain sev- 
eral well-separated jumps (bottommost loop in Fig.[l]^a)) 
indicating that the reversal is non-uniform and sequen- 
tial either because it becomes dominated by pinning of 
domain walls on defects or, equally likely, for intrinsic mi- 
cromagnetic reasons such as the bifurcations discussed in 
Ref. Il4 For \0h\> 88°, finally, sharp discontinuities are 
no longer visible in the experimental curves. Hereafter, 
we limit our discussion to angles less than 85°. 

The zero-current astroid of the free-element is shown 
in Fig. [U^b). It is qualitatively similar to the astroid of 
an isolated single-domain particle. However, two differ- 
ences may be pointed out. i) First, the experimental 
astroid is comparatively flatter: The switching field close 
to the hard-axis is markedly larger than that along the 
easy-axis. This is an indication that the lateral dimen- 
sions of the free-element are too large for truly coherent 
rotation of the magnetization to occur throughout the 
entire range of angles explored. Indeed, in the litera- 
ture, demonstrations of coherent rotation behavior exist 
only for smaller magnetic objects.'* Full micromagnetic 
simulations reproduce this flattening and show that, in 
the parallel state, an 5'-like non-uniform magnetization 
pattern forms in the free-layer at large field angles, ii) 
Second, unlike the square astroid of a Stoner particle, the 
astroid of the free-element shows a significant asymmetry 
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FIG. 2: Experimentally determined astroids of the free-element in a 100 nm x 200 rmi [Co/Ni]/Cu/[Co/Ni]/[Co/Pt] nanopillar 
for varying values of the DC current density: (a) j = +2 x lO" A/m^ (b) +3 x lO" A/m^, (c) -1-4 x lO" A/m^, (d) 
+6 X 10" A/m^ and (e), -6 x lO" A/m^ The zero-current astroid is also shown for comparison (solid grey symbols). 



with respect to the hard-axis of magnetization. This is 
due to the sizable dipolar coupling that exists between 
the two magnetic elements of the pillar only separated by 
a 4 nm thick copper spacer. The reference-element mag- 
netized along z produces an average stray field of about 
25 mT which tends to maintain the magnetization of the 
free-element oriented upwards. This results in a displace- 
ment of the astroid towards negative Hz , quite as for an 
exchange-biased ferromagnetic film.^^ 

We now extend these studies to include the contri- 
bution of spin-polarized currents to reversal. To avoid 
damage of the nanopillar during measurements, we have 
had to limit our investigations to moderate current den- 
sities (|j| < 6 x 10^^ A/m^). Moreover, the system 
studied behaves in such a way— that positive currents 
which become spin-polarized by transmission through 
the reference-layer are much more efficient at modify- 
ing the switching fields of the free-element (especially 
Hap^p) than negative ones polarized on reflection from 
the reference-layer. Therefore, we are mainly presenting 
here results obtained with positive currents. 

Positive currents favor an alignment of the free-element 
magnetization parallel to that of the reference-element. 
As may be seen in Fig. [2lja-d), the spin-transfer effect 
they generate modifies strongly the upper half of the as- 
troid, that is, the part of the critical curve correspond- 
ing to the switching from the antiparallel state to the 
parallel state [Hap^p)- As spin-transfer starts to op- 
erate, the easy-axis cusp present for zero DC current 
[Fig. [Ijb)] disappears. Instead, a dip forms along the 
easy-axis direction which gets deeper and broader as 
the current increases. At first, the large angle parts 
of the Hap^p branch of the astroid remains relatively 
unchanged [Fig. [Ua,b)]. At larger currents, the entire 
branch which can be probed with radial fields (con- 
stant Oh) is affected and takes a semi-circular shape 
[Fig. [2lc,d)]. In contrast, the lower half of the astroid 
corresponding to the reversal from the parallel state to 
the antiparallel state (Hp^ap) is not strongly affected 
under positive current. The domain of stability of the 
parallel state only slightly expands in the direction of 
negative Hz and the easy-axis cusp remains visible. 



For negative currents [Fig. [21[e)], a qualitatively sym- 
metrical behavior is observed. On the side of the Hp^ap 
branch, the astroid shrinks back as a dip forms along 
the easy-axis direction. On the side of the Hap^p 
branch, the astroid noticeably expands towards posi- 
tive Hz. Quantitative differences exist though (e.g., the 
size of the dip) between astroids obtained for currents 
of the same magnitude but opposite polarities (compare 
Figs-^d) and[2l[e)). These are straightforwardly related 
to the difference in the spin-transfer efficiency between 
the two current directions. The fact that upon injection 
of currents of the two polarities, large parts of the critical 
curve, and sometimes one half of it, remain virtually un- 
changed is a strong indication that the heat and Oersted 
field generated by passing the current through the device 
are not sufficient to modify the astroid and cannot be 
evoked to account for the observed distortionsjii^ 



In the remainder of the paper, we will mostly concen- 
trate on those parts of the astroid where spin-transfer 
gives rise to a reduction of the switching field, that is, in 
those regions where the dip forms, which are the most 
relevant to technological applications. Theoretical in- 
vestigations^? have shown that in the on- axis geometry 
{9h — 0), linear relations should exist between the val- 
ues of Hp^AP and Hap^p and the magnitude of the 
DC current injected, in partial agreement with experi- 
mental findings;^ The present results reveal that such is 
not the case if the field is applied at a large angle away 
from the easy-axis. More specifically, for every non-zero 
field angle, a threshold current jmin exists below which 
spin-transfer does not affect magnetization reversal. This 
minimum current increases with increasing 9h- As a con- 
sequence for moderate current values, there exist a field 
angle 0^"^ above which spin-transfer becomes ineffective. 
These results are rather counterintuitive. Indeed, as Oh 
increases, so does the relative angle between the mag- 
netization vectors of the reference and free layers and, 
consequently, so does the spin-transfer torque. Naively, 
one might therefore expect an enhanced efficiency of spin- 
transfer at large field angle. 
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FIG. 3: Calculated astroids of the free-element. The values of the DC current density used in the simulations are the same 
as experimentally (Fig.[2]|: (a) j ^ +2 x lO" A/m^ (b) -h3 x lO" A/m^ (c) +4 x 10^^ A/m^ (d) +& x lO" A/m^ and (e) 
—6 x 10^^ K/ra? . The large black symbols indicate the portions of the astroid accessible upon sweeping the magnetic field at 
constant Oh- The smaller grey symbols are the complementary parts determined from field sweeps at constant Hy. The red 
(respectively blue) lines delineate regions of the field space where the stationary antiparallel state (respectively parallel state) 
is replaced by a steady precession state with large negative (respectively positive) m^. The astroid computed for j = is also 
shown for comparison (grey line). 



III. NUMERICAL SIMULATIONS 

To gain a deeper understanding of the effect of spin- 
transfer on the astroid, we carried out numerical simu- 
lations in a macro-spin approach where the magnetiza- 
tion of the reference layer is supposed to be fixed. By 
comparing the experimental results with those from such 
a model, our goal is to unravel what is intrinsic to the 
physics of spin-transfer and what is possibly due to sam- 
ple imperfections, higher-order magnetic anisotropics or 
deviations from a uniform magnetization distribution. 

In order to include a description of the dipolar inter- 
actions in the system as realistic as possible under the 
assumption of uniform magnetization, the two magnetic 
elements of the nanopillar were assumed to be identical, 
200 nm long, 100 nm wide, 3 nm thick parallelepipeds, 
separated vertically by 4 nm. Analytical results from 
Newell et al}^ were used to calculate two tensors; first, 
the self demagnetizing tensor Nd , which relates the de- 
magnetization field Hf inside parallelepiped i to its mag- 
netization Mi through Hf = -Nd ■ Mi 
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0.9350 



(1) 



and second, the so-called mutual demagnetizing tensor 
Nm which allows one to give a simple expression for the 
stray field H^'^ generated by element j (magnetization 
Mj) and experienced by element i, H^'^ = —Nm ■ Mj 



Nm = 



0.0105 







0.0217 







-0.0322 



(2) 



Both the demagnetizing field and the stray field from 
the reference-element H'^'p were included in the effec- 
tive field acting the magnetization of the free-element. 



For the latter element, we assumed a saturation magne- 
tization of Ms = 650 kA/m, a perpendicular-to-plane 
magnetic anisotropy constant K± = 3 x 10^ J/m'^ and a 
damping parameter a — 0.01, whereas for the reference- 
element we used Mg"-'' = 500 kA/m (Ref. [H). 

The switching fields were extracted from magnetiza- 
tion versus field loops calculated by solving numerically, 
with the fourth-order Runge-Kutta algorithm, a modi- 
fied Landau-Lifschitz-Gilbert (LLG) equation including 
an additional spin-transfer torque acting on the magne- 
tization M of the form proposed by Slonczewski— li^, that 
is, —(3(9) j /Ms [M x (M x p)], with the direction of spin 
polarization p = z. The variation of the spin-transfer 
efficiency function (3 with the relative angle 9 between p 
and the magnetization direction m — M/Afs assumed in 
these calculations is of the same form as that found in 
Ref. [23 



/3(m-p) 



7?i 



1+ 



q- 



2dfj,oMse [60 + fci (m • p) 60 - &i (rh • p) 

(3) 

where d is the thickness of the free-layer, 7 > is the 
gyromagnetic ratio, and e > is the absolute value 
of the electron charge. The parameters q+ = 1/60, 
g_ = —1/600, 60 = 1, and bi — 1/2 were chosen in 
the following way. bo was arbitrarily set to 1. The ratio 
60 /bi was set to 2 so as to obtain a significant asymmetry 
of the slope of the spin-transfer-torque between angles 9 
close to 0° and those close to 180°. The ratio q+/q- was 
set to —10 to account for the asymmetry of the GMR 
device architecture. Finally, the ratio 60/9+ = 60 was 
adjusted so that values of the current density used in 
the simulations would approximately match those in the 
experiments. 

No effort was made to adjust these parameters further 
in order to get the best possible agreement with experi- 
mental data and reach a quantitative agreement, should 
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FIG. 4: Calculated change in the switching field for the transi- 
tion from the AP state to the P state, (a) AHap^pU, Oh) ~ 
Hap^pUjOh) — Hap^p{0,Oh) versus current density j for 
various values of the applied field angle 9h- (b) AHap^p 
versus field angle 9h for equally spaced values of the current 
density j ranging from -1-0.6 x 10^^ A/m^ (topmost, orange 
curve) to -f3.6 x 10^^ A/m^ (bottommost, black curve). 



this be possible. Yet, as may be seen in Fig. [31 most of 
the important features of the spin-transfer-distorted as- 
troid discussed before are accounted for by this approach. 
In particular, the existence of a threshold current when 
6h ^ is confirmed. This is even more clearly seen in 
Fig. lUa) which plots the change in the switching field 
AHap-^pU, 9h) = Hap^pU, 9h) - Hap^p{0, Oh) as a 
function of the current density, for various field orienta- 
tions. Also obvious from Fig. [H^a) is that the larger Oh 
the less linear the variation of the switching field with j, 
beyond jmin- Figure |DJb) illustrates the loss of efficiency 
of spin-transfer beyond an angle 6*^°^ which increases 
with increasing j . 

The absence of a large bubble-shaped extension of the 
astroid along the If^-axis, for negative current, is the 
only important point of disagreement between modeling 
[Fig- Ele)] and experiments [Fig. [^Je)]. Numerically, the 
extent of this profusion is found to be strongly reduced 
if the easy-axis of magnetization makes an angle of a few 
degrees with respect to the magnetic field plane. A slight 
distribution in the orientation of the easy-axis in the free 
layer of the studied device is thus a possible explanation 
for the discrepancy. We note however that the choice of 
the spin transfer efficiency function [3 affects also the size 
of this bubble. 



FIG. 5: (a) Polar and azimuthal angles 6 and </>, and vector 
basis {m, e^, e^}. 



IV. ANALYTICAL MODELING 

Interestingly, a simple extension of the Stoner- 
Wohlfarth (SW) model ^ is sufficient to capture most fea- 
tures of the current-distorted astroid. Let us consider a 
uniformly magnetized nano-magnet with saturation mag- 
netization Ms and uniaxial magnetic anisotropy of axis 
z and constant K {Hk = "^^K/ ^qAIs). Its magnetization 
M = Ms rii is described by a polar angle 9 measured from 
z and an azimuthal angle measured from x [Fig. [S] . As 
is the case experimentally, let us assume that the exter- 
nal field is applied in the yz plane. In the orthonormal 
direct basis formed by the three vectors m, eg — dm/ 89, 
^4> — {1/ sin9)drn/d(t>, the modified LLG equation that 
governs the dynamics of rii may be written as 



dm 
~(it 



-7(m X H*^^) + a(m X — ) (4) 



where a is Gilbert's damping constant. The total effec- 
tive field acting on M 



H:^.^ = H + HKirii ■ z)i+j ^^"^'^^ (m x p) 

is the sum of the usual effective field 

1 dE 



H 



^J.oMs dm. 

which derives from the magnetic energy 

E = -/ioM ■ H - K{m ■ if 
and of the spin-torque field 

HsT = j — (lii X z) = —j — sin 9 ba 

7 7 



(5) 



(6) 



(7) 



(8) 



At equilibrium, m is necessarily parallel to H*^^. 
Therefore, the equilibrium conditions are 

(H://.e,)„ = (9a) 
(H:^^-ea = (9b) 

where the subscript "0" denotes equilibrium. 
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The consequence of introducing spin-transfer is 
twofold, (i) First, since H^ff have a component 

along 60 [Eqs. I5|8j . rh does not always lie in the plane 
defined by the applied field H and the easy-axis z, at 
equilibrium. Therefore, the problem becomes three- 
dimensional, in general. In the limit of small currents 
(and/or small field angles), however, deviations of m 
from the (z, H) plane remain small and, in first order 
approximation, the problem can still be treated as if it 
were two-dimensional ((/) = 7r/2). (ii) Second, the sta- 
bility of the equilibria can no longer be determined from 
free energy considerations only and a new stability crite- 
rion must be derived. This is possible analytically in the 
small current limit. Indeed, performing a linear stability 
analysis around an equilibrium position 9o, i.e. relating 
stability to the gradient of the total torque along eg, leads 
to the criterion 



< 



(10) 




H/H 

y K 



y K 



FIG. 6: Comparison of astroids drawn from the extended 
Stoner-Wohlfarth model, Eqs. [TT] and 1121 (thick lines) and 
deduced from macro-spin simulations (symbols) in the case 
of a small positive current, j = -1-3 x 10^ A/m^ (a) and a 
small negative current, j — —3 x lO' A/m^ (b). As in Fig. |3] 
the thin red lines delineate regions of the field space where 
the antiparallel state is replaced by a steady precession state 
with large negative m^. 



The above criterion is written so as to reveal the compe- 
tition between two terms related to the torques produced 
by the effective field He// and the spin-torque field H5T, 
respectively. More practically, equation [10] shows also 
that stability in the off-axis case {9h 7^ 0, 6'o 7^ 0[7r]) is 
influenced not only by the spin-transfer efficiency func- 
tion I3{9) but also by its derivative 8(3/ 89, a point first 
put forward by Smith et al^ 

By combining the equilibrium condition, Eq. I9al and 
the stability criterion, Eq. 1101 we can derive the following 
parametric equations of the astroid under small current 



Hy = Hk s\t? 9o - j sin 9q C{9o) 
H^ = - Hk cos^ 6*0 - j cos9o C{9o) 



with 



Ci9o) 



1 8i(3sm9) 



89 



(11) 



(12) 



where the correction terms due to spin-transfer appear 
most clearly. As may be seen in Fig. [6][a,b), predictions 
from this analytical model regarding the stability of the 
equilibria closely agree with results from computer simu- 
lations using the same assumptions. Of course, the model 
is unable to predict the existence of dynamic states and 
regions of the field space where two magnetic states ex- 
ist but only one is stationary (as indicated by numerical 
simulations) are not comprised in the astroid thus de- 
fined. In the extended SW model, the negative counter- 
part of the dip is a kind of protruding bubble attached to 
the rest of the astroid through a crunode (double point, 
Fig.lSb)). 

For the numerical application of Eq. [11] and the com- 
parative macro-spin simulations shown in Fig. [6] the 
overall magnetic anisotropy constant K was set to — 
HqMs/I, with the values of and M5 given in sec- 
tion [ITT] Moreover, we used an expression of 13 in agree- 



ment with Slonczewski's original proposition- 



/3(m-p)- 



7?i 



2d^oMse 



1 



with 



fiP) 



/(P)(3-f m.p)-4 



(1 + n 



4p3/2 



(13) 



(14) 



and a degree of spin polarisation of the electrons coming 
out of the reference layer arbitrarily set to P = 50 %. For 
demonstration purposes, this function was deliberately 
chosen different from the one used in the numerical sim- 
ulations of Fig. [3] [Eq. [3] although both are monotonous, 
increasing functions of 9 on the [0,7r] interval. Yet, the 
two functions which were previously considered in the 
literature lead to qualitatively identical modifications of 
the astroid shape and to the existence of an angle 6*^5°^ 
above which the switching field is virtually unchanged. 
This demonstrates once more the robustness of this fea- 
ture. 



V. DISCUSSION 

The magnitude of the effective field acting on the mag- 
netization at the very beginning of reversal is the key 
parameter to consider in order to understand the effect 
of spin-transfer on the astroid. The reason for that is 
that Heff is a measure of the intrinsic damping. Indeed, 
the larger H^ff, the faster the precession of M around 
He//, in the event of an excursion away from the equilib- 
rium, and the stronger the viscous damping which drives 
M back to equilibrium. From Eqs. [Ha] and [TOl one can 
readily derive that along the zero-current astroid 



H, 



Hk sin 



(15) 
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The on-axis geometry is extremely particular. In this 
geometry, the effective field and the intrinsic damping 
vanish at the zero-current switching fields [Eq. [15] . This 
is why, close to these field values, any small amount of 
extrinsic damping brought in by spin-transfer affects the 
magnetization reversal. Furthermore, the orientation of 
the magnetization being independent of the applied field 
{9q = or tt), Heff varies linearly with H and the ex- 
trinsic damping due to spin-transfer is just proportional 
to j. This leads to the observed linear relation between 
the switching field and the injected current [Fig. Ufa)]. 

In the off-axis geometries, the situation is qualita- 
tively different. Indeed, the effective field always re- 
tains a sizable magnitude and the magnetization expe- 
riences a finite intrinsic (positive) damping at reversal. 
To induce an early switching of the magnetization, e.g. 
Hap^pU) < Hap^p{0), spin transfer must produce 
enough negative damping to overcome this finite intrin- 
sic damping, jmin is the smallest current density which 
realizes this. If j is less than jmin then, irrespective of 
j, reversal occurs when the equilibrium loses its local 
stability, i.e. upon crossing (exiting) the zero-current as- 
troid, as for j = 0. The mathematical complexity of the 
Oh ^ Q case, which is largely due to the fact that the ori- 
entation of the magnetization changes continuously with 
H , makes it difficult to derive an analytical expression 
for the dependence ol jmin on Oh- To explain why jmin 
increases with Or or, equivalently, why a moderate cur- 
rent density generates enough negative damping only up 
to a given field angle 0^"^ one has to invoke the fact that 
the intrinsic (positive) damping increases faster with Oh 
than the extrinsic (negative) one. 

VI. CONCLUSION 

In summary, we have investigated the effect of spin- 
transfer on the Stoner-Wohlfarth astroid of a small mag- 



net with uniaxial magnetic anisotropy. The distortions 
observed experimentally are well accounted for qualita- 
tively by both macro-spin numerical simulations and a 
simple extension of the SW analytical model. Evidence 
has been given that spin-transfer is more efficient at mod- 
ifying the switching field in geometries close to the so- 
called axial geometry— where the external field is ap- 
plied along the easy-axis and spin polarisation direction. 
On departing form this situation, a threshold current ap- 
pears below which spin-transfer is ineffective, the larger 
the field angle the larger this current. 

Our results have implications for solid-state device ap- 
plications. For zero-current, the SW model predicts that 
a two-fold reduction of the switching field can be achieved 
by applying the external field at an angle of 45° from the 
anisotropy axis, a strategy often used in magnetic mem- 
ory cells. What our results suggest is that combining such 
a strategy with spin-transfer will not necessarily help in 
reducing the switching field further. For the injection of 
current to be efficient, its density will have to exceed the 
threshold value jmin- We note finally that the physics 
which describes the distortions of the Stoner-Wohlfarth 
astroid under the influence of spin-polarized currents is 
similar to that which explains spin-transfer-induced in- 
stabilities of the sensing layer in OMR read heads of hard 
disk drives 
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